Swapping spin currents: interchanging spin and flow directions 
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We introduce a previously unknown spin-related transport phenomenon, consisting in a trans- 
formation [swapping) of spin currents, in which the spin direction and the direction of flow are 
interchanged. Swapping is due to the spin-orbit interaction in scattering. It originates from the 
correlation between the spin rotation and the scattering angle. This effect is more robust than the 
skew scattering, since it exists already in the first Born approximation. Swapping may lead to the 
spin accumulation with spin polarization perpendicular to the surface, unlike what happens in the 
spin Hall effect. 

"After all this is again a problem which falls into the same 
iry of problems related to the orientation of spin." 

From the report of our Referee E 



Spin-orbit interaction leads to the mutual transforma- 
tions of spin and charge currents [Il-Q : the charge current 
generates the transverse spin current and vice versa. The 
family of these phenomena includes the anomalous Hall 
effect in ferromagnets 0-01 and the direct and inverse 
spm Hall effect N, Q ■ After the first experimental obser- 
vations 0, Q , the spin Hall effect in semiconductors and 
metals attracted much interest and became the subject of 
numerous experimental and theoretical studies (stj.Most 
of the theoretical work is devoted to the intrinsic (not 
related to electron scattering) mechanism of spin-charge 
coupling joUlOll , which was first proposed by Karplus and 
Luttinger [6|. 

In this Letter we introduce another transport phe- 
nomenon belonging to the same family. We will show 
that because of spin-orbit interaction, a given spin cur- 
rent induces not only the charge current, but also a trans- 
verse spin current with interchanged spin direction and 
the direction of flow. We will refer to this transformation 
as the swapping of spin currents. 

The notion of spin current was introduced in Ref. [ij. 
It is described by a tensor qij where the first index in- 
dicates the direction of flow and the second one shows 
which component of spin is flowing. Following [s, 11[ we 



write down the phenomenological equations describing 
the coupling between spin and charge currents, qij and 
qi (more accurately, q is the electron flow density, related 
to the electric current density J by q = — j/e, where e is 
the elementary charge). We consider an isotropic mate- 
rial with inversion symmetry. Then we have [ill: 



FIG. 1: Schematics of electron scattering by a negative 
charge. The electron spin sees a magnetic field B ^ v x E 
perpendicular to the trajectory plane. Note that the magnetic 
field (and hence the sense of the spin rotation) has opposite 
directions for electrons scattered to the left and to the right. 



Pure symmetry considerations allow for additional 
terms in Eq. ([2|) proportional to and 5ijq^kk ^ which 
describe transformations of spin currents. In the presence 
of electric field E and spin polarization P, this would 
result in additional contributions to qij proportional to 
EjPi and 6ij{E ■ P), which were mentioned already in 

[Ell- 

However the physical origin of these contributions 
was not understood at the time and the effect was not 
studied ever since. 

We will show that the additional terms should always 



enter in a combination gj^'' —Sijq)^^! so that Eq. ^ should 
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be modified as: 
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where qf^^ and qf^^ are the primary currents, which may 
exist in the absence of spin-orbit interaction, Sijk is the 
unit antisymmetric tensor and 7 is a dimensionless pa- 
rameter proportional to the strength of spin-orbit inter- 
action. 



with a new dimensionless parameter yc (I4I . We will also 
show that the resulting swapping of spin currents origi- 
nates from the correlation between the scattering direc- 
tion and spin rotation during collisions. This effect is 
more robust than the spin-charge coupling: the swap- 
ping constant x exists already in the Born approxima- 
tion, while 7 appears only beyond this approximation. 
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The physical origin of the swapping effect can be read- 
ily seen from Fig. 1 illustrating the spin dependent scat- 
tering. The most important element is the magnetic field 
B existing in the electron's moving frame and seen by the 
electron spin. This field is perpendicular to the plane of 
the electron trajectory and has opposite signs for elec- 
trons moving to the right and to the left of the charged 
center. The Zeeman energy of the electron spin in this 
field is, in fact, the spin-orbit interaction. 

Three spin-dependent effects can be seen: 

1) The precession of the electron spin around B during 
the collision, leading to the Elliott- Yafet spin relaxation. 

2) The spin asymmetry in scattering (the Mott effect, 
or skew scattering) resulting from the additional force 
proportional to the gradient of electron Zeeman energy. 

3) The correlation between the directions of electron 
spin precession and scattering. Indeed, one can see that 
while the spin on the left trajectory is rotated clock- 
wise, the spin on the right trajectory is rotated coun- 
terclockwise. This correlation and its consequences for 
spin transport were not discussed previously. We will 
now see that this correlation leads to a transformation of 
spin currents. 

Suppose that the incoming electrons move in the y- 
direction and are polarized along y (spin current qlj^J). 
The electrons scattered to the left (right) will acquire 
a small positive (negative) projection of spin on the x 
axis. This means that the initial qy^J spin current is 
partly transformed to —qxx- For the case when incoming 
(along y) electrons are polarized along x, a similar rea- 
soning shows that the initial spin current q^} will give 
rise to qxy Thus in the latter case the spin direction 
and the direction of flow are interchanged. Summariz- 
ing, one can derive the spin current transformation law: 
(Ij^i^ ~ ^ijl'kk lij^ expressed by Eq. ([3]). 

The reason why the terms q^^^ and Sijq^i^ enter in such 
a combination can be explained as follows. Since the 
spin-orbit interaction is weak, the transformed currents 
are small compared to the primary currents. Suppose 
that the primary current contains only one component, 
say q^xx ■ In the first order in spin-orbit interaction there 
should be no correction to this component. (Similarly, 
when a vector a is rotated on a small angle, the first 
order correction Sa is perpendicular to a.) The combi- 



nation q^^^ — ^ij<tk^ is the only one to meet this require- 



ment. Generally, the terms q'^^ and —Sijq)^l! enter with 
coefficients whose difference is second order in spin-orbit 
interaction. 

We now proceed with the formal derivation of Eqs. ^ 
[3]). Introduce the spin density matrix Papip), where a 
and P are spin indices, p is the electron momentum. For a 
given energy, the evolution of pap due to elastic collisions 



,(0) 



is described by the following equation 0, [isj 



dpaf}{p) 

dt 



(4) 



where the integration is done over the angles of p' and 

'Z'{fP,^)*- (5) 



Here A'' is the impurity concentration, v = pjm is the 
electron velocity, and F'^^, is the scattering amplitude 
relating the initial state (a'p') and the final state (ap): 



pap 
a'p' 



(6) 



1? is a scattering angle, cr is the Pauli matrix vector and 
n is the unit vector perpendicular to the scattering plane 
n = p' X p/\p' X p\. 

The second term in the Eq. originates from the 
spin-orbit interaction Hso = A(fc x VU) ■ <t, where k — 
p/h, U{r) is the scattering potential, and A is the spin- 
orbit constant. 

In the absence of spin-orbit interaction B{{}) = and 
W^,% = Nv\A{§)\^6aa>Sf:)p>. Then Eq.© reduces to the 
conventional Boltzmann equation with the differential 
cross-section given by 1^(19)1^. 

Equation Q is applicable when the orbital motion can 
be considered classically It was previously used to 
study depolarization during atomic collisions [15| and the 
spin-charge coupling for J = 3/2 holes in the valence 
band as well as for the carriers in a gapless semiconductor 
[III, see also 

It is convenient to present p^^ in the form 



Papip) = T; [fiP)^aP + Tip) ■ CTafi\ , 



(7) 



where /(p) = Tr(/5) and T'(p) = Tr((T/5) are the particle 
and spin polarization distributions respectively. 

From Eq. (|4]) we obtain a system of coupled kinetic 
equations for /(p) and T'(p): 



rf/(p) _ 

dt 
dVjp) 
dt 



Nv dn' M/(p')-/(p))+a2n.p(p')], (8) 



= Nvjdn' [aiV{p') - T{p)) + n^nfip') 

+ asnxV{p')+<Tiny.{nxV{p'))], (9) 



where we have introduced four scattering cross-sections 
ai(i9) - + a2{^) - 2Re(AB*), 



CT3(i9) = 2Im(AB*), 0-4(1?) = 2|B| 



(10) 



As we will see, they have the following physical mean- 
ing. The usual transport cross-section is determined by 
CTi, containing a correction |-Bp, which is second order 
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in spin-orbit interaction. The spin-charge couphng, de- 
scribed by the parameter 7 in Eqs. dUIS]), is related to 0-2 
(skew scattering cross-section). The cross-section 0-3 is 
responsible for the swapping of spin currents. Finally, a4 
determines the Elliott- Yafet spin-relaxation rate. These 
cross-sections were, in fact, introduced by Mott and 
Massey [3] to describe the depolarization of an electron 
beam in a single scattering event. 

We note that in the Born approximation the ampli- 
tudes A{'d) and B{i}) have a phase difference n/2. Thus 
in this approximation (T2 =0, however the swapping 
cross-section (T3 is non-zero. 

Because of the rotational invariance of Eqs. ([51121), the 
system of these integral equations can be exactly solved 
for arbitrary initial distributions /o(p) and 'Po(p) [Hj 
16[. We will derive the coupled equations for the first 
moments of / and P defined as 



Q^-J^vJ{p), Q,^^J^v,V,ip). (11 



The charge and spin currents are given by 

q^^J d^p Q,(p), q^j = J d^p Qy (p). 



(12) 



We multiply Eqs. hy Vi and integrate over the 

angles of p. The integrals over dil in the right-hand sides 
can be easily evaluated. The result reads ^9|] 



(13) 



dQ^ 

Here the momentum relaxation time Tp — (Nvatr)^^ is 
expressed in the usual way through the transport cross- 
section (Tfr 



dn cri(l?)(l - COS!?). 



(15) 



The parameters 7, >f, and the Elliott- Yafet spin relax- 
ation time Ts are given by 



1 



7 = 



2(Ttr 
1 

2atr 



drj cr2(i?)sini?, (16) 
dnasi^) sin^, (17) 



Nv / dn cr4(l?). 



(18) 



Introducing the averaged over angles of p spin polariza- 
tion P{s), in a similar way we obtain 



dPje) _ P{e) 
dt T, 



(19) 



For a given energy e = p^/(2to), Eqs. [HI de- 
scribe the decay of the charge and spin currents and the 
spin polarization due to elastic collisions. 



We are interested in the stationary solutions of Eqs. 
(|131ll4p in the presence of external driving forces, for ex- 
ample in the presence of an electric field. The right-hand 
sides of these equations will acquire additional terms Q^f^ 
and Qij"*, which are proportional to the derivative dfo/ds, 
where /o(e) is the equilibrium distribution function. For 
degenerate electrons, dfo/de oc S{e — Ep), where Ep is 
the Fermi energy. Integrating the stationary solutions of 
Eqs. (fT51 [TH) over p we finally derive Eqs. (jl] where 
the parameters 7 and >f are taken at the Fermi energy. 
We note that Eqs. (HI O (but not Eqs. dUEl)) are 
valid only in the first order in spin-orbit interaction [l^] ■ 
For non-degenerate electrons, the expression (|T7l) should 
be modified as 



N 



(73 sin {) dCl 



(20) 



where (...) means averaging over the electron energy with 
the equilibrium distribution function. 

We will now calculate the swapping parameter >c for 
the Coulomb scattering in the Born approximation. For 
a positively charged center the scattering amplitudes are 



(21) 



where as = h^e/me^ is the Bohr radius, e, m and e 
are the dielectric constant, electron effective mass, and 
charge respectively. Using Eqs. ([TUl [T71) . and (HH) we 
obtain 



X = 2Xk^ 



(22) 



In semiconductors with the band structure of GaAs, 
in the limit of small effective mass, the Kane model 
gives [iOl A = h^/AmEg for A > E^g and A = 
{h'^/3mEg){A/Eg) for A < Eg, where Eg is the band 
gap and A is the spin-orbit splitting of the valence band. 
Hence 



Ef 

>f = — for A > Eg, 
^9 

'iEp A ^ 
K = for A < Eg. 

3 Eg Eg 



(23) 
(24) 



For a bulk electron concentration of lO^^cm ^ we cal- 
culate a quite large value >c = 0.3 for InSb (Eq. ([23l) ) 
and X = 0.003 for GaAs (Eq. ([241) 1. 

To observe the swapping effect one should obviously 
have a primary spin current. Most easily this is achieved 
by passing an electric current through a ferromagnet. 
However the strong exchange magnetic field will tend 
to destroy spin currents with polarization different from 
that of the majority spins. Whether the swapping effect 
in ferromagnets can produce a measurable effect, or not, 
is an interesting question, which we do not address here. 
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(b) Primary spin current Secondary spin currents 




FIG. 2: Schematics of a proposed experiment revealing the 
swapping of spin currents, a - the qyy spin current is electri- 
cally injected in a semiconductor (S) through a ferromagnetic 
contact (F) with a magnetization along the y axis. The wavy 
arrow symbolizes optical detection of the 2;— component of 
spin polarization near the surface, b - the swapping effect 
transforms the primary spin current qf^J into Qxx and q^^. 
This should lead to the appearance of excess spin polariza- 
tion at the lateral boundaries of the sample (Pz < at the 
top face and > on the bottom one). This polarization 
may be detected optically. 



A possible way to see the swapping effect in semicon- 
ductors, where the exchange field is negligible, is pre- 
sented in Fig. [2I The primary spin current qlj^J is pro- 
duced by spin injection in a semiconductor 3D sample 
through a ferromagnetic contact. Swapping will result in 
the appearance of transverse spin currents Qxx = Qzz — 
—xq^y. Those secondary currents will lead to an ex- 
cess spin polarization near the lateral boundaries of the 
semiconductor sample, which could be detected by op- 
tical means. At the top face there will be a polariza- 
tion Pz < Q and at the bottom face Pz > 0. (Similarly 
Pj; < at the front face and P^; > at the back face.) 
The accumulation of spins polarized perpendicular to the 
surface distinguishes this manifestation of swapping from 
the spin Hall effect, where the accumulated spins are par- 
allel to the surface. 

It seems plausible that not only the scattering mech- 
anism considered here, but any mechanism leading to 
the spin Hall effect, will also give rize to the swapping 
phenomenon. However for the moment we do not see 
how this could happen due to the "intrinsic" mechanism 
This problem should be addressed in future work. 

In summary, we have considered a new phenomenon 
caused by spin-orbit interaction and consisting in swap- 



ping of spin currents. This effect originates from the 
correlation between the spin rotation and the direction 
of scattering. 
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the Russian Ministry of Education and Science, and the 
President Grant for young scientists (MD-1717.2009.2). 
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APPENDIX. REPORT OF REFEREE B 

The authors claim to have found a new "spin current 
swapping effect" based mostly on symmetry considera- 
tions and obtaining a result of polarization at the edges 
due to a polarized charge current in linear order to the 
spin-orbit coupling. 

I find the result problematic from several considera- 
tions (see below), with some of the explanations and 
some of the partial conclusions already present in dif- 
ferent forms in other publications (although the authors 
seems to be, according to the authors, the only active 
spintronic researchers in the field), and therefore cannot 
recommend the manuscript for publication in Physical 
Review or Physical Review Letters. 

The authors, although not specified in the abstract, 
restrict their study to the case of weak spin-orbit cou- 
pling (otherwise the notion of spin-current and their sym- 
metries considerations do no longer hold) and identify 
several tensor components in the equations relating spin 
current and charge current and identify a relating term 
Kappa which they try to calculate. 

They make a qualitative argument for the effect based 
on the previously used explanations of the work of Sinova 
ct al 04 (not cited here) which identify the rotation from 
effects in between collisions. Note also that similar cou- 
plings as the ones mentioned have been obtained in the 
linear response theory of other authors (e.g. Sinitsyn et 



al PRB 081212 (04)); also, many of their results I believe 
to be nothing more than the Edelstein effect which has 
the same symmetry (but again, not cited by the authors). 

The authors follow the procedures of Ref. 11 (which is 
mostly equivalent to the result of Ref. 1 , simply reprinted 
30 years later in Physical Review Letters) and seem to 
derive (with many assumptions along the way) Eq. 20 
which they obtain the final result Eq. 22 which they pro- 
ceed to estimate for GaAs and InSb. One obvious prob- 
lem is that their result has not been observed the experi- 
ments (ref. 7) and according to their geometry it should 
have been observed but of course it was not, suggest- 
ing strongly that in fact x = 0. Such spin-accumulation 
would have been observed, very easily according to their 
estimates, in these experiments. A linear result in SO 
coupling in these case should be highly suspicious since 
it should have been extremely evident a long time ago. 

I also find it intriguing that the authors suggest that 
experiments 7 and 8 generated a large interest and debate 
when in fact in 2004 a very large debate was ignited by 
prior theory works, one of which the authors use similar 
diagrams (Fig. 1) in the heuristic explanations of their 
argued effect. 

Given all the inconsistencies that I find in their ar- 
gument I cannot agree at this point to recommend the 
submitted manuscript to Physical Review and much less 
to Physical Review Letters. 



